Entanglement in helium 
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Using a configuration-interaction variational method, we accurately compute the reduced, single- 
electron von Neumann entropy for several low-energy, singlet and triplet eigenstates of helium atom. 
We estimate the amount of electron-electron orbital entanglement for such eigenstates and show that 
it decays with energy. 
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I. INTRODUCTION 

The role of entanglement as a resource in quantum 
communication and computation [l], 0] has stimulated 
many studies trying to unveil its fundamental aspects 
as well as to provide quantitative entanglement mea- 
sures [3-5]. More recently, the role of entanglement at- 
tracts growing interest in systems relevant for chemistry 
and biology, for instance in light-harvesting complexes, 
governing the photosynthesis mechanism in proteins [6|- 
\T\\ . In such instances, due to the complexity of the in- 
volved systems, ab initio treatments are computationally 
prohibitive and one is forced to consider simplified (spin) 
models, with parameters introduced phenomenologically. 

Besides computational problems for many-body quan- 
tum systems, one must address the problem of measur- 
ing entanglement for indistinguishable particles fl2l - [l6j ]. 
The main difficulty in quantifying entanglement is due 
to the symmetrization or antisymmetrization of the wave 
function for bosons or fermions. That is, to discriminate 
entanglement from correlations simply due to statistics 
of indistinguishable particles. In spite of this difficulty, 
bipartite entanglement has been investigated in a num- 
ber of systems of physical interest, including the entan- 
glement dynamics of electron-electron scattering in low- 
dimensional semiconductor systems fl7| and the entan- 
glement of low-energy eigenstates of helium atom [l8| . 

To evaluate entanglement in helium, Ref. [IH used 
high-quality, state-of-the art Kinoshita-type wavefunc- 
expressed in terms of Hylleraas coordi- 



and then computed the purity of the 
one-electron density operator by means of a 



tions 
nates 
reduced 

twelve-dimensional Monte Carlo numerical integration. 
In the present paper, we compute several low-energy 
states of helium by means of a simple configuration- 
interaction variational method. Our approach has sev- 
eral advantages. First of all, we do not need to evaluate 
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multidimensional integrals: the reduced density matrix is 
obtained by purely algebraic methods. The reduced den- 
sity matrix can then be easily diagonalized and therefore 
we can access not only the linear entropy but also the 
von Neumann entropy. Finally, we express our varia- 
tional, Slater-type basis, in terms of (radial and angular) 
single-particle coordinates, and therefore such Fock-state 
basis could be easily extended to many-body systems. 
Despite the above mentioned, still unsolved conceptual 
difficulties in the definition of entanglement for indistin- 
guishable particles, we propose a way to evaluate the or- 
bital entanglement for states close to Fock states. From 
such a measure we conclude that the amount of entan- 
glement exhibited by helium eigenstates | Is, ns; 1 S) and 
| Is, ns; 3 5) drops with energy. 

The paper is organized as follows. In Sec. |TT] we de- 
fine the measure used to estimate the entanglement of 
helium eigenstates. In Sec. IIIII we shortly describe the 
variational method that we used to determine the low- 
energy reduced density operator p\. Our numerical re- 
sults for entanglement of the helium atom are discussed 
in Sec. IIVI Finally, our conclusions are drawn in Sec. [V] 



II. ENTANGLEMENT OF HELIUM 
EIGENSTATES 

The non-relativistic Hamiltonian of the helium atom 
reads, in atomic units, 
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where Z = 2 denotes the nuclear charge, Pi the mo- 
mentum of electron i {i = 1,2), r,; its distance from the 
nucleus and r 12 the inter-electronic separation. 

Since we are neglecting the spin-orbit interaction, we 
can consider global wavefunctions S factorized into the 
product of a coordinate wavefunction ^Sf and a spin wave- 
function x : 



H CTliCT2 (ri,r 2 ) = *(ri,r 2 )x ( r ic 



(2) 
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The overall state must be antisymmetric and therefore 
a measure 12j of the amount of entanglement E(\E)) 



of the state 3 has been proposed in terms of the von 
Neumann entropy of the reduced density operator R\ — 
Tr2(|S)(S|) of one patrticle, say particle 1, obtained after 
tracing the overall, two-body density operator over the 
other particle: 



£(|S)) = S{Rx) - 1, 
with the von Neumann entropy 

S(iii) = -^A i log 2 A i 



(3) 



(4) 



where {A^} are the eigenvalues of R\. 

However, with such definition a first problem arises. 
When considering the triplet subspace, spanned by the 
spin states xtt> 7f (Xn+Xtf), and Xii> it is clear that we 
should consider the case -^=(xti + Xit) separately from 
the cases xtt and Xii- Indeed, -^(xti + Xit) is a maxi- 
mally entangled Bell state of the two spins while xtt and 
Xn are separable states. Therefore, the standard spec- 
troscopic characterization in terms of triplet and singlet 
states is no longer useful for the purposes of computing 
entanglement and one should study separately the entan- 
glement properties of the states composing the triplet. In 
this context, we would like to point out that, neglecting 
spin-spin interaction, the choice of the basis states span- 
ning the triplet state is completely arbitrary and that the 
above discussed {xtf, 75 (xtl + Xit), Xii} is onl y one in 
between the infinite possible choices. 

To avoid this ambiguity, in this paper we compute the 
entanglement for the orbital part * only of the wavefunc- 
tion. Since * is antisymmetric for spins in the triplet 
subspace, we can measure the amount of entanglement 
of * as follows: 



where 



EQ*)) = 5(pi) - 1, 



S(Pl) = -^Z Ai l0g 2 Al 



(5) 



(6) 



is the von Neumann entropy of the reduced density op- 
erator pi = Tr2(|*)(*|), and {Aj} are the eigenvalues of 
Pi- 

When the spin part of the wavefunction is in the sin- 
glet state xs = 75 (xti - Xit), the orbital part is nec- 
essarily symmetric and this causes an additional, open 
problem in the quantification of entanglement. Indeed in 
this case the reduced von Neumann entropy alone is not 
sufficient to discriminate between entangled and separa- 
ble states [15] . The core of the problem is the fact that 
we can have separable states with either S(pi) = or 
S(pi) = 1. The first instance corresponds to basis states 
of the kind ^u(ri, r 2 ) = 0i(ri)(/>i(r 2 ), the second to basis 
states like #y(ri,r 2 ) = ^[<^(ri)^-(r 2 ) + 0j(ri)0;(r 2 )], 



with i 7^ j. On the other hand, even the quantification of 
entanglement of the global, antisymmetric wavefunction 
by means of Eq. ((3]) poses a problem. Indeed, as the von 
Neummann entropy is additive for tensor products, for 
the state 5 = ^uxs Eq. © gives E(\E)) = 1, while for 
the state S = &ijXS (* 7^ j) we have E(\E)) = 2. Even 
though measure ([3]) gives different results, the amount 
of entanglement in both cases is the same, since the or- 
bital wavefunctions and are both separable: the 
reduced density matrices for these two states have differ- 
ent entropies only due to the symmetrization of the state 

We will not address in this paper the unsolved prob- 
lem of quantification of entanglement for a generic state 
*. On the other hand, since from our calculations it 
turns out that for each low-energy helium eigenstate the 
reduced density operator for the orbital part is rather 
weakly perturbed with respect to one of the two above 
non-entangled cases, we expect that an approximate 
quantification of entanglement is provided by the dis- 
tance between the von Neumann entropy S{p\) of p\ 



and the von Neumann entropy S(p^) (S(pi'') = or 
S(pi) = 1) of the reduced density operator p± for the 



corresponding non-interacting, non-entangled state: 
E(\*)) = \S(p 1 )-S(p? ) )\. 



(7) 



Note that for antisymmetric orbital wave functions this 
definition reduces to ([5]). 



III. METHOD 

We compute with high accuracy the lowest energy 
states of helium by means of a variational method, 
the configuration-interaction method (see, for instance, 
Ref. [HI)- Orthonormal basis functions are provided by 



,h ,m 1 ;n 2 ,l 2 m 2 ( r l > r 2) 

F ni j i;n2t i 2 (ri,r2)Yi imi (0 1 )y; 2m2 (f2 2 ), 



(8) 



where Y/ m are spherical harmonics, with solid angle for 
the particle i and the radial functions £ , n 1 ,/ 1 ;n 2 .b( r ii r 2) 
are obtained after orthonormalizing the Slater-type or- 
bitals 



Rai{r) =r n+ '- 1 exp(-e„, i r), 



(9) 



with £ n> i variational parameters, and properly sym- 
metrizing the products of the obtained one-particle radial 
wafefunctions f n i(j*). That is, if the spin wavefunction is 
symmetric, F must be antisymmetric, 



-fn2h(n)fnih(r2)]; 



(10) 
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if the spin wavefunction is antisymmetric, F must be 
symmetric: 

Fn 1 ,l ll n2,h( r l,r 2 ) = -J=[.f ni h( r l)fn 2 l 2 (r2) 

v2 (11) 

+fn 2 i 2 {ri)f ni i 1 {r 2 )] 
if (ni,/i) ^ (n 2 ,l 2 ), 

F ni ,h;ni,h( r l> r 2) = /ni2i( r l)/Wi(»"2) (12) 

otherwise. 

We then compute the reduced (single-electron) density 
matrix 

tt(ri,ri)= /dra^ri.raJ^Cri.ra), (13) 

with 

*(ri,r 2 ) = ^ c /lj2 $ /li/2 (n,r 2 ), (14) 

/l,/2 

with the multi-indexes I\ = (ni,li,mi) and 1% = 
{n 2 ,l 2 ,m 2 ). Since the expansion is done over an or- 
thonormal basis the reduced density matrix on that basis 
is simply given by a partial trace over the second particle 
of the overall density matrix: (pi)^,/; = J2i 2 PhJrJ'^hi 
where pi 1 j 2 -j> i j^ = ci 1 i 2 c* i , r . We point out a major 
advantage of the configuration-interaction method: the 
reduced density matrix is obtained by purely algebraic 
methods, without numerical computation of the multi- 
dimensional integrals of Eq. (fl"3|). The reduced density 
matrix can then be easily diagonalized and in this paper 
we will study the entanglement properties of helium by 
means of the eigenvalues {A&} of p\. 

IV. NUMERICAL RESULTS 

In the following, we will present data only for total or- 
bital momentum quantum number L = 0, thus implying 
lx = lv = l and mi = — m 2 . 

We first discuss convergence of our method, as a func- 
tion of the number n max (?) of radial wavefunctions for a 
given I and as a function of the cut-off Z max on I. For 
the low-energy states discussed below, we found that 
! max = 2 (S, P, and D shells) and n max m 10 — 20 (from 
"max = 10 for the ground state up to ri max = 20 for the 
higher excited states reported below) are sufficient to re- 
produce helium eigenergies with at least four significant 
digits (as deduced from comparison of our results with 
those of Refs. [H dH-il) and reduced von Neumann 
entropy S'(pi), estimating at least two-three significant 
digits. To illustrate the convergence of our method, we 
provide in Table U the obtained values of S(pi) of the 
ground state of helium for different values of the cut-offs 
( max and n max (we take the same n max for all values of 

0- 



rimax = 5 0.04131 0.07772 0.07844 0.07833 

n max = 6 0.04133 0.07776 0.07848 0.07837 

n max = 10 0.04134 0.07777 0.07849 0.07839 

rimax = 11 0.04134 0.07777 0.07849 0.07839 



TABLE I. Reduced von Neumann entropy of the ground state 
of helium, computed with different cut-off values in the basis 
of Slater-type orbitals. 



State 


S( P i) 


Si (pi) 


l(Ls)VS) 


0.0785 


0.01606 


lls^s; 1 ^) 


0.991099 


0.48871 


lls^s; 1 ^) 


0.998513 


0.49724 


lls^s; 1 ^) 


0.999577 


0.49892 


lls^s; 1 ^) 


0.999838 


0.499465 


lls^s; 1 ^) 


0.999923 


0.499665 


lls^s; 1 ^) 


0.999961 


0.499777 



TABLE II. Reduced von Neumann and linear entropies for 
the lowest energy singlet eigenstates of helium. 

The reduced von Neumann entropy S(pi), as well as 
the linearized entropy Sl(pi) = 1 — Tr(pi) often used in 
the literature, are shown in Table iLTl and in Table IPLT1 for 
several low-energy singlet and triplet eigenstates, respec- 
tively. 

Since the obtained values of the von Neumann entropy 
are very close to those expected for Fock states, which 
are separable, the entanglement content is weak and can 
be estimated by means of Eq. ([7]). The obtained results 
are shown in Fig. Q] as a function of the state number 
n, for both singlet states \ls,ns; S) and triplet states 
| Is, ns; 3 S). Note that data, with the exception of the 
ground state value of entanglement, are consistent with a 
power law decay of entanglement with n. From a power- 
law fit we obtained E(n) — 0.19 n~ 4 - 41 for singlet states 
at n > 2 and E(n) = 0.040 n -319 for triplet states. The 
same entanglement data are shown as a function of en- 
ergy in Fig. [2] It can be clearly seen that the entangle- 
ment content drops with energy. This result is rather in- 
tuitive in that for states {lSjUS'^S} and |ls, rts; 3 S') the 
wave functions corresponding to the states Is and ns 
are localized far apart for large n. Therefore, electron- 



State 


S(pi) 


Si (pi) 


\ls,2s; 3 S) 


1.00494 


0.500378 


|ls,3s; 3 S> 


1.00114 


0.5000736 


lls^s; 3 ^} 


1.000453 


0.5000267 


|ls,5s; 3 S> 


1.000229 


0.5000127 


|ls,6s; 3 5> 


1.000133 


0.5000070 


lls^s; 3 ^} 


1.000091 


0.5000047 



TABLE III. Same as in Table [TT1 but for the lowest energy 
triplet eigenstates of helium. 
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electron interactions become weaker (and entanglement 
smaller) when n increases. We note that, due to a dif- 
ferent definition of entanglement, the energy-dependence 
obtained in our results contrasts the one obtained in 
Ref. [Hj]. 




FIG. 1. Entanglement of the singlet states | Is, ns; 1 S) (circles) 
and of the triplet states | Is, ns; 3 S) (triangles). 




w 



FIG. 2. Same as in Fig. [TJ but with the entanglement of the 
helium eigenstates plotted as a function of their energy W 
(measured in hartrees). 



V. CONCLUSIONS 

In this paper, we have computed the reduced von Neu- 
mann entropy for several low-energy (singlet and triplet) 
eigenstates of helium. The von Neumann entropy has 
then been used to estimate the amount of entanglement 
of such states, showing that the entanglement of the 
states |ls, ns; ^ and llSjns; 3 ^) decays with n, that is, 
with energy. This result is in agreement with the intu- 
ition, suggesting that when the electronic wavefunctions 
are localized far apart, electron-electron interactions are 
weak and therefore the entanglement is expected to be 
small. Our results can be seen as one of the first steps 
towards a "spectroscopy of entanglement" for atomic and 
molecular systems. While in this quest helium atom 
constitutes the simplest example, the variational scheme 
used in this paper may be extended to more complex sys- 
tems, and the obtained numerical results could be used 
to test the validity of perturbative approaches or of sim- 
plified phenomenological models. Finally, we point out 
that, while our bipartite entanglement measure can be 
applied to wavefunctions close to separable Fock states, 
quantification of entanglement for generic states of indis- 
tinguishable particles remains an interesting open prob- 
lem. 
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